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Abstract 

We study the electroweak phase transition by lattice simulations of 
an effective 3-dimensional theory, for a Higgs mass of about 70 GeV. 
Exploiting, among others, a variant of the equal weight criterion of phase 
equilibrium, we obtain transition temperature, latent heat and surface 
tension, and compare with Mh ~ 35 GeV. In the broken phase masses 
and Higgs condensates are compared to perturbation theory. For the 
symmetric phase, bound state masses and the static force are determined. 
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1 Introduction 



Recent lattice studies of the electroweak phase transition had been trig- 

gered by the interest in understanding the baryon number asymmetry of the uni- 
verse. Significant cosmological consequences would require a sufficiently strong 
first order transition ||, |f(J. The transition has to be strong enough, both in 
order to accomplish a sufficient rate of baryon generation during the transi- 
tion and to prevent the wash-out of baryon number after it is completed. The 
present quantitative understanding of possible mechanisms as well as experi- 
mental lower bounds for the Higgs mass make this unlikely within the minimal 
standard model. Extensions, in particular supersymmetric ones may still be 
viable, however Q-(f3|. 

A second reason for lattice investigations was the wish to control the be- 
haviour of perturbative calculations of the effective action. This quantity is 
the appropriate tool of (non-lattice) thermal quantum field theory for dealing 
with symmetry breaking. Infrared problems prevent a perturbative evaluation 
of the free energy in the symmetric phase to higher loop order. However, the 
true, non-perturbative nature of the symmetric phase will be characterized by 
massive W— and Higgs bound states instead of massless W gauge bosons. A 
self-consistent approach to provide masses across the transition, e.g. by gap 
equations ]l4| , can improve the ability to calculate perturbatively the symmet- 
ric phase. Gauge field condensates are another property of the symmetric phase, 
expected to lower its free energy density. 

An independent method is needed to characterize the electroweak phase 
transition even within a pure SU(2) gauge-Higgs version of the theory. This 
model has become a test-field to control the validity of perturbative predictions 
over a broad range of Higgs masses. At present, one is interested to see whether 
the first order transition ends somewhere around a Higgs mass Mh ~ 100 GeV 
[[l5j fUSf , Lattice simulations are not only able to describe both phases starting 
from first principles but, moreover, make it possible to put both phases into 
coexistence near the phase equilibrium. Thus one is able to measure directly 
quantities like latent heat, surface tension, condensates etc. quantifying the 
strength of the transition. 

One approach to lattice calculations of the electroweak transition is based 
on an effective 3-dimensional Higgs model. It is attractive phenomcnologically 
because it circumvents the problem of putting chiral fermions on the lattice. 
Due to dimensional reduction, fermions as well as non-static bosonic modes 
contribute to the effective action. In contrast to QCD, dimensional reduction 
should work for the electroweak theory around and above the transition tem- 
perature because g 2 is small. For the electroweak phase transition this approach 
has been pioneered by Farakos et al. (see e.g. |l7|, |l^|). This program aims 
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at exploring the accuracy of dimensional reduction at various Higgs masses by 
comparing various parameters of the transition with those of 4-dimensional lat- 
tice and perturbative approaches. Perturbation theory is necessary to relate 
the 4-dimensional continuum theory to the parameters of the dimensionally re- 
duced theory and, finally, to the bare coupling parameters of the lattice action. 
Dimensionally reduced versions retain the remnant of the temporal gauge field 
A (as an adjoint Higgs field) or not (as in this work). 

In its simplest version the dimensionally reduced effective theory is again a 
SU(2)-Higgs theory with just one doublet. This model represents whole classes 
of id theories but might not be sufficient to describe all Ad variants generically. 
So far, only this simple effective theory has been studied by lattice techniques, 
for several Higgs self-couplings A3 in units of the 3d gauge couplings squared 
[0- H If- Here we extend our analysis Q to a higher coupling value 

A3 1, M* 2 

g\ 8 l 80GcV J [) 



namely A3/P3 ~ 0.095703, referred to as M'^ = 70 GeV, and compare with 
previous results at this smaller coupling Xz/g^ ~ 0.023926 (M^ — 35 GeV). As 
expected the first order nature has become weaker but is still evident. 



2 Relation of the dimensionally reduced model 
to the SU(2) Higgs theory in four dimensions 

The procedure of dimensional reduction in gauge theories consists of two steps. 
In a first step, nonstatic modes of the gauge field are integrated while the tem- 
poral component A$ is retained as adjoint Higgs field. The action is further 
simplified by integrating over Aq in a second step. The mass of the latter field 
is given by the Debye mass tod. The reduced Higgs theory has the following 
action resembling the form of the 4-dimensional theory 




where a, (3,b = 1, 2, 3. The 3-dimensional mass parameter 7713(^x3), renormalized 
at scale ^3 and the renormalization group invariant couplings A3 and gs can be 
expressed by dimensional mapping in terms of the running parameters of the 
4-dimensional theory at the scale fj, and temperature T in the MS scheme. The 
renormalization scale \i is most conveniently chosen as fi — fix = 7.05T which 
puts logarithmic contributions of bosonic origin to zero. 
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To make this paper self-contained we collect some formulae for later use. 
As the result of the second step of reduction mentioned above, the couplings 
<73 and A3 and the 3-dimensional mass 7713 can be expressed in terms of the 
corresponding couplings and mass of the intermediate theory with Aq included, 
g 3 , A3 and 7713 (and the Debye mass mj = ^/5/6gT). The latter parameters 
in their turn can be written in terms of the temperature and the parameters of 
the 4-dimensional theory: 

2 

2 _2/i 93 \ 

93 = 33U 



247T77J£) 



g\ = g 2 ^ T )T(l + 9 



16tt 2 3' 



T 3 ^ 

A3 — A3 — 



1287TTO 



D 



n 2 SM 2 

j a _ A (w) r(i + I f-^), (3) 



■> 3 _ 2 g 3 /15 /x 3 

m|(M3) = - 2 (^) + T(iA 3 + ^) + -Z!_(^ g4 + ^ 3 2 ) 
' 8 ^ + 9A333 " 12A^)(log — - 0.348725). 



167r 2V 16 yd ™ A M/v b M3 



The 4-dimensional renormalized quantities g 2 {n) , ^ 2 (m) and A(/z) have the 
generic form with one-loop corrections formally indicated 



A(/i T ) 
v 2 {Ht) 



Sg 2 ) 



9o 2 (l 

^(1 

8 M 2 / 



5\) 



(4) 



with g 2 = 4:V2G F M^. G f = 1.16639 x 10~ 5 Gey~ 2 is Fermi's constant. The 
corrections indicated in these formulae, which depend on g 2 ,/iT and the Higgs 
and W mass squared Mfj and M^, can be found in Ref. ttq]. For definiteness, 
we use as the 3-dimensional renormalization scale /i3 =53 and take Mw = 
Mz = 80.6 GeV. For the 4-dimensional coupling in the loop corrections we 
choose g 2 = g 2 . 
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At Mjj — 70 GeV the corrections in ([!]) without fermions are numerically 
obtained as follows 

5g 2 « -0.01343- 0.01945 log(|u|/M^), 
5X w -0.01625 + 0.008491 log(/i|/M^), 

(5^ 2 w -0.01892- 0.004721 log(/i|/M^). (5) 

Similarly, we find for M* H = 35 GeV 

5g 2 w -0.01281- 0.01945 log(M|/M^), 
<5A w -0.05396 + 0.06131 log(/i|/M^), 

^ 2 w -0.07286- 0.005766 log(/i|/M^). (6) 



3 The 3— dimensional lattice model 



3.1 Lattice action 

On the lattice, we study the St/(2)-Higgs system with one complex Higgs dou- 
blet of variable modulus. The gauge field is represented by unitary 2x2 link 
matrices U x , a an d the Higgs fields are written as $> x = p x V x . p x = ^Tr($+$ x ) 
is the Higgs modulus squared, V x an element of the group SU(2), U p denotes 
the SU(2) plaquette matrix. The lattice action is 



S = f3 G J2(l-lTrU p )-p H Y,lTr(<PtU x , a <Z> x+a ) 
p I 

+ -I) 3 ) (7) 



(summed over plaquettes p, links I and sites x), with 

/fc = A- (8) 
The lattice Higgs self-coupling is 



91 PG 
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and the hopping parameter can be expressed in the form 



2(1 - 20h) 
4 



= tl 2 2 - (10) 

6 + a z mi 



The relation between the 3-dimensional bare mass squared m|, expressed via 
the lattice couplings (see (|lo|)) and the renormalized continuum mass squared 
has been worked out by Laine [jl9| to two loops 

to 3(a*3) = m\ + m\+m\ 

m\ = ^(^53+6A 3 ), E = 3.175911 (11) 

ml = I ^((^ 3 3 4 + 9A3. 9 3 2 -12A2)(log-^- + 0.09) 
+5.0^ + 5 .2A 3ff3 2 ). 

This completes the relation of the lattice parameters of the 3d model to the Ad 
continuum parameters. 

The lattice model defined by (Q) is numerically studied at given couplings 
Pg ) Ph and A3 / . We have used the same simulation algorithms as in our pre- 
vious investigations. We have combined a 3-dimensional Gaussian heat bath 
update for the gauge fields with a 4-dimensional Gaussian heat bath method 
for the Higgs field. In the last one, the Gaussian step was improved for later 
acceptance by taking the non-linearity into account. To reduce the autocorre- 
lations near to the phase transition, a heat bath step as described was followed 
by several reflections (eight in practice) for the Higgs (and partly one reflection 
for the gauge field). 

In the search for the phase transition, bulk variables like the space-averaged 
square of the Higgs modulus and the space-average link 



2 x \ 2 

P = T T T >,Px 



Ei 



lilt k 



- T 

•v L z x 

3mE^W, (12) 
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are heavily used (the average quartic Higgs modulus p is defined analogously) . 
Others like the average plaquette 



P = 



3L x L y L z 



G 



are less useful for this purpose, although they also show a two-state signal on 
large enough lattices. In principle, the histograms and discontinuities of all 
bulk variables have to be known in order to estimate the strength of the phase 
transition, but p 2 proves to be the most important one. It has turned out that, 
due to the confinement property of the symmetric phase, histograms of medium 
size rectangular Wilson loops show a two-state signal as well. 

While the bulk quantities discussed are measured after each Monte Carlo 
iteration, other quantities like Wilson loops and correlation functions are mea- 
sured after every 10th iteration. 

To reduce the noise in the evaluation of Wilson loops the original links are 
replaced by their mean field values ||(J. They are calculated from the links 
interacting with the original links through the staples. In the case of SU(2) this 
value is calculated analytically. This procedure is not applied to the links in the 
corners of the Wilson loops. 

The Higgs and vector boson masses are measured from correlation functions 
of extended operators of length n (r b are the Pauli matrices) 

1" : ^„W = \Tr(T b ^+U x ^JJ x+a . a ...U x+{n _ l)a ^ x+na ) (14) 

For the mass fits we have used the maximum n(= 4) for the vector boson which 
gave the best signal in the symmetric phase. In the Higgs channel we did not 
observe such a strong n dependence. In the analysis we used both n = and 
n = 4. 

3.2 Implementation on QUADRICS parallel computers 

The simulations which are reported in this paper have been performed on 
QUADRICS parallel computers at the University of Bielefeld, at the QUADRICS 
Q4o and the CRAY-YMP of HLRZ in Jiilich. Additionally, codes have been 
partly tested on the Ql at DESY-IfH Zeuthen. The computing facilities at 
Bielefeld were provided by the Deutsche Forschungsgemeinschaft (DFG) to the 
groups being part of the DFG task force program " Dynamical Fermions" . The 
code written in TAO had to be portable to various topologies of the QUADRICS 
family (Q4o, QH2) of parallel machines. The Q4o is a 32 node machine with a 
2x4x4 topology, the QH2 has 256 nodes in 4 x 8 x 8 topology. 

For the calculation of large Wilson loops we have defined long strings of 
links and their mean field values (as auxiliary field attached to the point of 
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beginning). The factors must be communicated from remote processor nodes. 
Forming loops in a later step, these auxiliary building blocks must be multiplied 
non-locally where nodes far away have to be addressed simultaneously. Due to 
memory constraints (size of the auxiliary field) we were able to measure Wilson 
loops on the Q4o for lattice sizes up to 40 3 . Similar communications have to be 
considered for correlation functions on non-neighbouring nodes. 

The accumulation of histograms of bulk variables as well as Wilson loops 
must be done on the host computer (after taking the appropriate global sums 
which are then available on all nodes simultaneously) since integer arithmetics 
is not possible on the QUADRICS processors. 



3.3 Multihistogram technique and phase separation 

As usual, the search for the phase transition point requires extensive application 
of the multihistogram technique [lj] to the bulk variables listed above. We 
have processed data obtained in up to 7 runs with different (3h values (for a 
given volume) within and near to the metastable region. Each run consisted 
out of 30000 to 145000 measurements, depending on the measured integrated 
autocorrelation time Ti nt . The respective maximum of the autocorrelation time 
increases roughly linearly with the volume. These numbers, observed in the run 
nearest to the respective pseudocritical /3# c , range from 200 (for 30 3 ) to 1500 
(for 64 3 ). These are autocorrelation times with tunneling and mainly indicate 
the longer lifetime of each one of the two metastable phases, i.e. the suppression 
of tunneling in the larger lattices. 

We look for the phase transition point in terms of the critical hopping pa- 
rameter, (3hc, for various (in fact two, in the present work) values of the lattice 
gauge coupling f3c, while the continuum couplings are kept fixed. In other 
words, we have studied the phase transition driven by m^. Then the lattice 
Higgs self-coupling (3r varies with (3h (see @). Therefore, the reweighting uses 
not only En n k, but p 2 and p A at the same time. In general, a 3-dimensional 
binning has to be performed in the two relevant parts of the action (correspond- 
ing to (3h and and some other observable of interest. This enables to create 
histograms for any bulk variable and for any value of (3h near to the transition 
point, based on an estimated density of states which subsums all measurements 
in the metastable region with appropriate weights. 

We have determined the finite volume pseudocritical 0h c (L) by the minima 
of the Binder cumulants 

B Elink (L,p H ) = l- {Ehnki) (15) 
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and B p i[L 1 [3h), by the maxima of the susceptibilities 

C Elink (L,p H ) = (E lmk 2 ) - (E Hnk ) 2 (16) 

and C p i{L, Ph), and using the equal weight method. Our aim is to apply this 
method for phase transitions of relatively weak first order. For these transitions 
a frequent tunneling between the pure phases is observed in the metastability 
region for all lattice sizes studied. This has made necessary the refinements 
described below. 

Concerning histogram methods at first order transitions, there is an ambi- 
guity in the literature how to define the critical coupling (in fact: pseudocritical 
for finite volume) . Throughout the metastable region, a two-state signal is vis- 
ible in the histograms, e.g. of p 2 and Eu n k, but there are two prescriptions for 
pHc- equal height of the maxima vs. equal weight of the two competing phases. 
While theoretical arguments favour the equal weight criterion, this requires 
in practice a procedure to separate the (measured or reweighted) histogram into 
contributions from configurations to be attributed to the respective pure phases. 
In addition, there are inhomogeneous (mixed) configurations contributing to the 
histograms. There is no generally accepted procedure that unambiguously de- 
fines the weights of these three contributions. 

Our main assumption is that the pure phases can be described by Gaussian 
distributions for any volume averaged quantity. Order parameters like p 2 or 
Eunk can be considered in this context. We have utilized the equal weight 
method as described below for the case of p 2 . The normalized histogram has 
been presented as a weighted sum of three histograms referring to the two pure 
phases and to all inhomogeneous mixed states 

p{p 2 1 (i H ) = W b p b (p 2 ,/3 H ) +W s p s (p 2 ,f3 H ) +W m ixPmix(p 2 ,PH) (17) 

with 

Wb + W s + W m i x = 1. (18) 

Wb, s denotes the weight of the broken/symmetric phase, w m i X is the correspond- 
ing weight of the mixed state, all weights are Ph dependent. The pseudocritical 
Phc(L) is then found for wj, — w s < 0.5. 

The positions, widths and weights of the pure phase histograms at a given 
Ph have been obtained by fitting the outer flanks of the two-peak histogram 
to Gaussian shape. At the same time, this fixes the weight w m i X and the p 2 
distribution to be attributed to configurations with domains of both phases in 
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equilibrium. We should mention potential additional problems to this fit proce- 
dure when histograms of both phases are strongly overlapping. In asymmetric 
transitions (generic for the Higgs case), with susceptibilities (widths) very dif- 
ferent in both phases, this happens for smaller lattice sizes. For sufficiently large 
volumes one can avoid this problem. 

We have used an iterative procedure to find the critical (3hc according to the 
requirement Wb =w s . It consists of merging all data into a single histogram at 
a tentative (3h- This step is followed by a fitting procedure as described above 
which tells the weights of the pure phases at this 0h- If they are not equal, (5r 
is corrected accordingly. 

Phase separation of measurements is necessary not only for this application 
of the multihistogram method. Another particular example is the measurement 
of pure phase correlation lengths near to the phase transition. For this purpose 
and for the splitting of badly separating histograms in some particular variables 
it is better to use the Monte Carlo time sequence of configurations for runs 
in the metastability range. The aim is to remove successful tunneling escapes 
and unsuccessful tunneling attempts towards the "wrong" phase from what 
should then be considered as the Monte Carlo trajectory restricted to the "right" 
phase. The procedure rescans the records of the volume averaged p 2 which has 
a well separated two-peak signal for all considered volumes. Referring to this 
variable a lower cut for the upper (broken) phase and an upper cut for the lower 
(symmetric) phase can be chosen. These cuts are determined in such a way that 
the remaining histograms (for the "pure" phases) are almost symmetric around 
their maxima. If the Monte Carlo history of p 2 enters the range of a certain 
phase and stays there for more than 100 iterations, all measured quantities 
(recorded for the whole trajectory, including correlators etc.) are considered to 
contribute to the statistics of the given phase until the trajectory leaves this 
phase again. 

The minimal time the trajectory is required to stay within one phase has 
to be much larger than the autocorrelation time without tunneling but smaller 
than the autocorrelation time with tunneling. The choice of 100 iterations is 
consistent with this requirement. 



4 Localization of the phase transition 

In Tables |l| and || the statistics is reported for each set of couplings and volumes, 
that has been used in the various procedures for the localization of the phase 
transition. 
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30 a 


48 a 


64 a 


0.343440 




75000 




0.343480 


30000 






0.343520 




85000 




0.343540 


50000 


80000 




0.343544 




120000 




0.343546 






90000 


0.343548 




145000 


120000 


0.343550 






100000 


0.343560 




40000 




0.343580 




110000 




0.343600 


40000 







Table 1: Data samples used for the localization of the transition for = 70 
GeV and (3 G = 12 





32 a 


40 a 


48 a 


0.340700 






45000 


0.340780 


40000 


40000 


45000 


0.340800 


40000 


100000 


90000 


0.343820 


40000 


40000 


45000 



Table 2: Data samples used for the localization of the transition for = 70 
GeV and /3 G = 16 



4.1 The infinite volume limit for 0h c 

In order to determine the set of pseudocritical values of the couplings @h (and 
the corresponding /3r) we have used the volume averages p 2 and Eu n k searching 
for the minima of the Binder cumulants (|l5| ) B p 2(L, (3h) and BE link (L, 0h), and 



for the maxima of the susceptibilities (16) C p 2(L, (3h) and CE link (L, (3h)- In 



addition to these methods we have used the equal weight method in the variant 
described above. 

To demonstrate the two-peak structure we show in Fig. [l] the measured 
histogram of p 2 on lattices 48 3 and 64 3 , all at (3q = 12, for (3h values nearest 
to the respective pseudocritical (3h c (L). The positions of the maxima change 
already only slightly with the volume. 

In Figs. H and || results of the multihistogram interpolation of our data 
for 0g = 12 for the Binder cumulant of Eu n k and for the susceptibility of p 2 
are presented. Finiteness and shrinking of the Binder cumulant with increasing 
volume present evidence for the first order nature of the transition at Higgs 
mass Mft = 70 GeV. The maximum of the interpolated susceptibility (with the 
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Figure 2: Multihistogram interpolation of BE link , (3g = 12 



background susceptibility subtracted) increases almost linearly with the volume. 

A histogram reweighted to pseudocritical (3h c {L) = 0.3435441 (as deter- 
mined according to the equal weight criterion) for a lattice size 64 3 is presented 
in Fig. U together with the Gaussians describing the pure phases just at that 
13 He- The distribution attributed to mixed configurations with domains of both 
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0.3435 0.3436 



Figure 3: Multihistogram interpolation of C p i , @g = 12 



i i 




lattice with contributions of the pure phases 

phases in equilibrium is well identified between the two peaks. 

The various pseudocritical (3h c {L) values for the three methods applied to 
Eii n k and p 2 are collected in Figs. || and ||. They are plotted there versus l/L 3 
for Pg = 12 and 16. The errors of the individual pseudocritical (3h values 
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Figure 5: Infinite volume extrapolation of (3h c for /3g = 12 
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Figure 6: Infinite volume extrapolation of (3h c for @g = 16 

were obtained by a variant of the Jackknife procedure. From each of the Monte 
Carlo histories near to the critical point (n simulated (3h values, n = 3 in 
most cases) one half has been left out of the analysis. This amounts to 2™ 
Jackknife estimators for [3h c - The equal weight method could be reasonably 
applied only for the larger volumes because of the overlap problem mentioned 
above. Corresponding to each method, a 1/L 3 fit has been used to yield a 
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respective /?|f c . The extrapolations nicely coincide as expected. 

In table || the extrapolations for each method are collected together with 
the average /3% c for j3 G = 12 and 16. With the use of formulae (0,|n],|j|) the 
Phc 1S translated into a physical temperature and an "exact" Higgs mass Mh- 
For definiteness, these numbers are given for the case of the SU (2) Higgs theory 
without fermions. Comparing these temperatures there seems to be not much 
space left for 0(a) corrections. The "exact" Higgs mass is practically the same. 





B 


C 


w b = w s 


ft; = 12 




Elink 


0.3435434 


0.3435430 




P 2 


0.3435429 


0.3435429 


0.3435441 


Wo 




0.3435433(6) 




T c /GeV 




154.34(1) 




Mjj/GeV 




66.52 




ft; = 16 




Elink 


0.3407950 


0.3407937 




P 2 


0.3407943 


0.3407939 




Wc 




0.3407942(6) 




T c /GeV 




154.68(2) 




Mh 1 GeV 




66.52 




Table 3 


Infinite volume limit for (3h c at M H 


= 70 GeV 



For comparison, at the smaller coupling (M H = 35 GeV) the transition 
temperature is T c = 85.7(1) GeV with the Higgs mass M H = 33.0 GeV. This 
has been obtained for gauge couplings in the range from Pa = 12 to 20 on 
lattices of size 40 3 and 20 3 . 



5 The strength of the phase transition 

5.1 Condensate discontinuities at the phase transition 

The jumps in {p 2 } and (p 4 ) are connected to the renormalization group invariant 
discontinuities of the quadratic and quartic Higgs condensates. The two-state 
signal for (p 2 ) and (p 4 ) is still clearly visible for all lattice sizes considered at 
the higher Higgs mass of M H = 70 GeV, where the transition turns out much 
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weaker than at Mjj = 35 GeV. The continuum condensate jumps can be put 
into relation to the lattice quantities by the following formulae 

A^+^/gj = ^f3 G (3 Hc A(p 2 ) (19) 
A((0+0) 2 )/ ff 4 = (i/3 G /fe c ) 2 A(p 4 ). (20) 




Figure 7: Infinite volume extrapolation of A(p 2 ) 

Corresponding to the different criteria applied for the definition of the pseu- 
docritical /3h c (L) we obtain histograms of the various operators just at the 
respective pseudo-criticality. From the jump of the operator expectation values 
between the phases (read off from the maxima of the corresponding histogram) 
an infinite volume extrapolation is performed. A collection of discontinuities 
of (p 2 ) and (p 4 ) for various finite lattices is shown in Figs. ^ and ||. As in 
the analysis of Ref. |?J we found that the size dependence of the jumps for all 
available lattice sizes is best described by a 1/L 2 fit. The extrapolations to in- 
finite volume (for the two criteria B p 2 and C p i ) and the averages (of the results 
obtained with different methods to extrapolate) of the jumps of the scalar Higgs 
operators are given in table ^ in lattice units. 

The corresponding condensate discontinuities in continuum units are col- 
lected in table ||. The quadratic Higgs condensate is already independent of 
finite a effects. On the contrary, the quartic condensate shows a severe a depen- 
dence. So we conclude, it is more subtle to extract an appropriate continuum 
value for this higher condensate. 



l(i 




Figure 8: Infinite volume extrapolation of A(p ) 







Cp 2 


average 






/3g = 12 


0.490(9) 


0.479(8) 


0.485(6) 


Pg = 16 


0.372(8) 


0.362(7) 


0.367(5) 


A(p 4 > 




/9g = 12 


4.86(9) 


4.81(8) 


4.84(6) 


/fc = 16 


3.62(8) 


3.51(7) 


3.56(5) 



Table 4: Infinite volume limit for A(p 2 ) and A(p 4 ) at M* H = 70 GeV 







A((^0) 2 )/3 3 4 


ft? = 12 
A? = 16 


0.250(3) 
0.250(4) 


1.28(2) 
1.65(3) 



Table 5: The continuum Higgs condensate discontinuities at = 70 GeV 



For the case of the lighter Higgs mass M# = 35 GeV we present the cor- 
responding quantities in table ^[ These data were obtained from two separate, 
completely metastable runs on a 40 3 lattice at 0g = 12 and (3h — 0.34140. 
They failed to tunnel to the other phase. The actual (3h value of these two runs 
had been found as pseudocritical one on a smaller lattice. 

Besides of jumps in the quadratic and quartic Higgs condensates also a 
discontinuity in the expectation values of Eu n k appears being a good indicator 



17 



A(^) 
A(p 4 ) 


6.24(1) 
96.1(2) 


A<(0+^) 2 )/ 5 | 


3.20(1) 
25.2 (1) 



Table 6: A(p 2 } and A(p 4 ) and condensates at M* H = 35 GeV and f3 G = 12 



for the phase transition as well. 

It is known (see jl7| ) , that the sum of expectation values of the various Higgs 
operators is constant according to a sort of Schwinger-Dyson equation 

- 3/3 H (Eunk) + (1 - 2f3 R )(p 2 ) + 2/3 R (p 4 ) = C (21) 

independent of the couplings (3h, @g an d A3/53. We have checked this sum rule 
for all used coupling values and found C = 2 exactly within very good numerical 
accuracy. 

From eq. (|2l]) one sees that the Higgs condensate discontinuities are related 
to each other by the following sum rule 

- 3(3 H A(E lmk ) + (1 - 2/3 R )A(p 2 } + 2(3 R A(p 4 ) = 0. (22) 

Additionally, the expectation value of the average plaquette (P) shows a 
discontinuity as well at the phase transition. The jump in this observable is 
numerically a tiny effect at the larger Higgs mass (mostly superimposed by the 
fluctuations, see the discussion and Fig. |l0| below). Nevertheless, we are able to 
estimate this jump using the phase separation technique discussed earlier. 

In table the jump A(P) is reported for /3h values nearest to the critical 
ones at = 70 and 35 GeV at lattice sizes 64 3 and 40 3 , respectively. At the 





A(P) 


M* H = 70 GeV p G = 12 
ft; = 16 


0.00037 
0.00015 


Mjj = 35 GeV f3 G = 12 


0.00370 



Table 7: Estimated plaquette jump A(P) 



larger M* H the phase separation technique is used. We do not set these jumps 
in correspondence to a continuum gauge condensate discontinuity 

A{\F« F^)/gt = -3^| A(P) (23) 
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since severe a effects and perturbative contributions are expected. 



5.2 Latent heat 

The latent heat Lheat = Ae (e is the density of thermal energy) of the transition 
is calculated according to |l7| 

^ = ^A(cPU) (24) 



With the reported jumps and the values for Mh and T c (g 2 (nT c ) ~ 0.38) we 
find the latent heats at the higher and smaller Higgs masses as given in table |L 





Lheat /T^ 


M* H = 70 GeV 
Mjj = 35 GeV 


0.0176(3) 
0.180(1) 



Table 8: Latent heat 



The equal weight method makes it possible to reconstruct directly the free 
energy densities of the pure phases in the vicinity of the phase equilibrium. The 
corresponding numbers are obtained in the iterative search of the pseudocritical 
(3 He using the reweighting technique at fixed [3g- The latent heat can then be 
expressed alternatively as the jump Ae of the energy density by 



J heat 



1 



T c 4 T c 2 a 3 L 3 dT 
Using at fixed f3a 



( log w s - log w b 



T=T a 



(25) 



d_ 

dT 



T=T a 



a*M* H d 
fHc 2T C 8P H 



0h=Pe 



(26) 



one finds 

Lheat 



si I 

T3 8 



(27) 



The change of the weights very close to the critical (3h is shown in Fig. ||. 
Taking the corresponding slopes needed in d27]) we obtain at Mjj = 70 GeV 
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for the largest lattices (where the equal weight method can be applied without 
ambiguity) L heat /T? = 0.0180(9) which is in good agreement with the other 
method. 



-3.4-1 0" 6 
-3.5-10" 6 
-3.6-1 0" 6 
-3.7-1 0" 6 
-3.8-1 CT 6 

3.8- 1CI 6 4.0- 10" 6 4.2- 1CT 6 

P H - .3435400 

Figure 9: Weights for the pure phases around the equal weight pseudocritical 
P H for p G — 12 on a 64 3 lattice 



5.3 Wilson loops at the phase transition 

If the lattice is large enough, a two-state signal becomes visible also in the 
average plaquette P. This can be seen in Fig. [l(] for (3g — 12, (3h = 0.343548 
at 64 3 . 120.000 measurements are collected in this histogram. 

For medium size rectangular Wilson loops a two-state signal is even better 
visible than for the plaquette. This is due to the fact that the configuration 
averages follow an area and perimeter law in the symmetric and broken phase, 
respectively. If these loops are not too large (i.e. not too small numerically), 
this can be well separated in the histograms of Wilson loops as demonstrated in 
Fig. [ll]. In that figure a histogram of 2.000 measurements of a space-averaged 
Wilson loop of size 12 x 12 is shown obtained during 20.000 MC iterations in the 
metastability region at lattice size 48 3 . The corresponding average plaquette for 
this size does not show a two-state signal. 

In Fig. [l^ we plot the logarithm of expectation values of symmetric L x L 
Wilson loops versus area A — L 2 , separated into symmetric and broken phase 
contributions at the critical (3h, for both Higgs masses = 70 GeV and 35 
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Figure 10: Measured histograms of P for fie = 12 
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0.015 0.020 
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Figure 11: Measured histogram of the 12 x 12 Wilson loop 



GeV (in both cases for /3g — 12). The numbers of measurements are 2.000 in 
both cases. This figure shows clearly the area law in the symmetric phase and 
practically no (or only weak) dependence of the string tension (in lattice units) 
on the Higgs self-coupling A3. 
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O M M =35 GeV, symm. phase 
- M H =35 GeV, broken phase 
□ M H =70 GeV, symm. phase 



-15_ ■ M H =70 GeV, broken phase 



100 200 300 400 

A 



Figure 12: Wilson loops for (3q = 12 



5.4 Estimate of the surface tension 

The coexistence of both phases in lattice configurations opens the possibility 
to determine the surface tension a. The use of the equal weight method as 
described above allows to estimate the contribution of the mixed phase state at 
the pseudocritical coupling for which Wb = w s (for large enough lattices). The 
mixed phase state weight w m ix is directly related to a. Using these weights to 
estimate a seems to be more natural than to obtain the surface tension from 
the ratio of the maximum to the minimum of the p 2 distribution at equilibrium 
defined by equal weight. For asymmetric transitions the choice of the maximum 
is even somewhat arbitrary (unless equal height is the criterium of choice). 

To get an idea about the shape of the interfaces (bubbles, walls) we show a 
snapshot of a particular configuration on a 32 2 x 128 {Mjj — 70 GeV, $q = 12) 
in the pseudocritical (3h region. This configuration has a p 2 = 3.67 where 
mixed states can be expected. To suppress the large fluctuations in the local 
p 2 values we have averaged them iteratively over the next neighbours. The 
configuration obtained in this way is shown in Fig. [l^ where we have plotted 
all lattice points with p 2 > 3.67. The filled (empty) region can be interpreted 
as the broken (symmetric) phase contribution. Besides of a few bubbles we see 
two large interfaces separating the pure phases. Although the configuration has 
been smoothed the interfaces are rather structured which makes it difficult to 
assign them a definite area A. For simplicity we use in the fit below A = L 2 a 2 , 
which tends to overestimate a. 
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Figure 13: Smoothed mixed state configuration, plotted are the points to be 
assigned to the broken phase 



We parametrize the relation between the weights at pseudocriticality and a 
for lattices of the form L 2 , x L, as follows 



= b L\ \o%L x exp(-2aa 2 L 2 x /T c ). (28) 



With that ansatz we try to describe data obtained either from cubic and from 
cylindrical lattice geometries. 

The factor L\ is an entropy factor which accounts for the positions of the two 
surfaces, logL x is the result (for d — 3) of the capillary wave approximation for 
the (internal) fluctuations of the surfaces. The degeneracy factor b is introduced 
to count different possible orientations of the surfaces dominant for cubic (b = 3) 
and for prolongated lattices (b = 1). 

In Fig. |l4| we present the data for a/T c 3 as function of (flc/L 2 .) for various 
lattice sizes and geometries and different (3c values at Mjj = 70 GeV. From a 
linear fit to the infinite volume limit we find the upper bound 

^ pa 0.00022. (29) 

^ c 

These lattice results are smaller by one order of magnitude than the one- 
loop estimates for the surface tension |23)| . A trend of finding smaller surface 
tensions at larger Higgs masses was already observed in . Whereas the latent 
heat is determined by the position of the broken minimum alone, the surface 
tension is sensitive to the shape of the effective potential in the whole tp range 
between the symmetric and the broken phase. The disagreement between the 
measured surface tension and the one-loop estimate at intermediate Higgs mass 
appears to indicate that the loop expansion to the effective potential gets out 



23 




of control at intermediate ip values already, reflecting the infrared problems of 
the symmetric phase. 

The surface tension is the only quantity indicating a substantial deviation 
from perturbative predictions at Mjj = 70 GeV, i.e. additional weakening of 
the phase transition with increasing Higgs mass. At still larger Higgs masses 
{Mh > M\v = 80 GeV) the phase transition has been suggested recently to 
be of second order jl5) or to be a smooth crossover jl6| in accordance with 
predictions from the study of gap equations Jl4|] . The first case would be difficult 
to discriminate from a phase transition being very weakly first order. In any 
case, however, one will observe the turn-over to a very weak or continuous 
transition at higher Higgs mass. 



6 Broken phase and perturbation theory 

Analysing the lattice data one should answer the question whether the broken 
phase can be understood perturbatively. There are no infrared problems in this 
phase because the elementary excitations are massive. Nevertheless, a pertur- 
bative treatment is expected to break down at larger Higgs masses, close to the 
phase transition. 

It is known p3 that for calculating the temperature dependent effective 
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potential the appropriate effective expansion parameter for the broken phase is 



2 

93eff 



9l 



2mw(T) 



(30) 



with the 3-dimensional gauge boson mass m\y(T). This coupling can be deter- 
mined from the two-loop effective potential |25| . Its value depends on the gauge 
fixing as well as on the choice of the renormalization point ^3. Using Feynman 
gauge and /X3 = g 2 for Mjj — 70 GeV this effective coupling is found to be 1.10 
at the critical temperature. 

In Figs. [l5| and [l6| we present a comparison of our lattice data with two- 
loop continuum predictions for the renormalized Higgs condensate 2 (T) and the 
vector boson (mw(T)) and Higgs boson (mij(T)) masses at the corresponding 
temperature, respectively, using the effective potential ^f|. The temperature T 
is implicitly given via m 2 (g 2 ) (see eqs. (|[ ||| ||). The {(f 2 ) data are obtained 
from the measured (p 2 ) by subtracting the one-loop and two-loop counterterms 

m 
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Figure 15: Comparison of perturbative and lattice 



4? = 2(0+0)( 33 2 ), 
(0 + 0>(5 3 2 )/3 3 2 = £/fc/MpV§^-^(log|/fc + 0.67). (31) 

We differ from B by a factor of 2 in the convention of defining 4> 2 . The two- 
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Figure 16: Measured masses at MJj = 70 GeV and 0g = 12 compared to 
continuum perturbation theory 

loop predictions are calculated in Feynman gauge. In the case of mw and 
tuh predictions are shown with and without the one-loop wave function renor- 
malization constants for the Higgs and vector boson masses Zh and Zw p4| 
(m = Z _1 / 2 too)- The prediction for tuh is derived from the curvature of the 
effective two-loop potential at the broken minimum, which is not identical to 
the pole mass. However, the difference is expected to be small. 

Comparing perturbation theory predictions with the data we observe very 
good agreement deeper in the broken phase, as should be expected. Wave func- 
tion renormalization is obviously required. Otherwise the agreement is poor. 
Close to the phase transition we observe a systematic difference between lattice 
data and the continuum calculation as function of m§. Nevertheless, the pre- 
dicted mass and (</> 2 ) values at the corresponding continuum critical point agree 
well with the measured ones at T c (m| < 0). The mapping of the two-loop 
results from the parameter v /v used in |2|] to m\{g^) fails close to the phase 
transition. This may be an indication that higher loop terms start to play a 
significant role in this regime. 



7 Some properties of the symmetric phase 

3-dimensional SU(2) pure gauge theory is known to possess confinement, i.e. 
the spectrum is formed by massive VF-balls instead of the elementary massless 
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Ws. Adding a Higgs doublet has similar consequences as adding fermions: new 
massive (0 + </>) bound states occur, and the static confining potential should be 
screened. The Lagrangian mass of the scalar bosons increases with temperature, 
i.e. sufficiently away from the phase transition the Higgs bound states become 
heavy. As a consequence, the symmetric phase should more and more resemble 
pure SU(2) gauge theory at higher and higher temperature. In the (3h range 
(resp. 777.3 range) that we have explored this could not yet be confirmed, however. 
The lowest Higgs bound state is still significantly lighter than the lightest 
ball. 



7.1 Higgs and vector boson bound states 



Results for the lowest Higgs bound state (0 + ) and the vector boson bound state 
(1~) are shown in Fig. [I?] as function of m^(g^)/ g^. Data from different (3q 
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Figure 17: Masses in 
to p G = 16 



the symmetric phase for (3q = 12, diamonds correspond 



values nicely coincide. Our results for the vector boson mass rrt\y should be 
considered as upper bounds, because of a possible admixture of states with a 
somewhat higher mass. For the Higgs mass mn this problem is not severe 
because the mass gap to higher mass states is significantly larger. This problem 
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has been carefully studied recently by Philipsen et al. |2(|. Note that our 
gauge boson masses are calculated from correlators of gauge invariant operators 
in contrast to Ref. p5[ where gauge variant correlators have been analysed, 
obtained in conjunction with gauge hxing (Landau gauge). 

One important conclusion from Fig. |l^ is the scaling behaviour of masses 
when plotted versus m|. Outside the immediate vicinity of the phase transition 
no A3 dependence is seen within errors. Our data point at the lighter Higgs 
mass has been reanalysed as compared to our Ref. Q . 

An attempt has been made to reproduce these bound state masses by approx- 
imate analytical methods [p7j based on the Feynman-Schwinger representation 
of correlators. The measured static potential (see below) is used as input. The 
level spacings (2s— Is, lp—ls) are very well reproduced by this approach for the 
data point of Ref. |^6| . The other interesting phenomenon is the m| dependence 
of the bound state masses, also being reproduced well. 



7.2 The static potential 

In four dimensions one usually considers Creutz ratios of Wilson loop expecta- 
tion values instead of the loops directly in order to avoid explicit renormaliza- 
tion. This is not necessary in three dimensions because renormalization becomes 
very simple. All ultraviolet divergent contributions are covered by the exponen- 
tiated (Abelian) one-loop term, having a logarithmic divergence. For the poten- 
tial, the one-loop massless T^-exchange contribution is gf 3/(87r) log(i?/(2a)). 
This is the result of a continuum calculation with a introduced as ultraviolet 
cut-off. Interesting enough, the behaviour of the two-loop contribution can sim- 
ply be predicted on dimensional grounds to be C2ig\R, with some constant C21 
not yet evaluated. Correspondingly, higher loop orders generate higher powers 
of R. Perturbation theory will therefore be appropriate to describe small Wil- 
son loops, i.e. the potential at small R, whereas it breaks down in a power-like 
manner at large distances (infrared regime). This is different from the 4c? case 
where corresponding effects are logarithmic. 

At small R, an appropriate ht to the potential (in units of <?f) should be 
provided by 

V(R)=V 1 + ^-\og(^) + c 1 R, (32) 

where a is identified as the lattice constant, and V\ has been introduced as 
additional fit parameter. It could be fixed, however, by a one-loop lattice cal- 
culation (not yet done). The renormalized potential will anyway contain some 
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arbitrariness in the overall normalization, due to the choice of the renormaliza- 
tion point. 

A priori, it is not clear up to which distance this ansatz may be appropriate. 
Interesting enough we shall find (see below) that it describes our data reason- 
ably well in the whole range, up to about R = 6/gf. Two-loop perturbation 
theory generates a term looking like a string tension. The value of ci fitted 
at intermediate R may therefore have perturbative as well as non-perturbative 
contributions. Because of the gauge coupling being dimcnsionful it will be dif- 
ficult not only in practice but also conceptually to separate perturbative from 
non-perturbative effects. This in particular concerns the definition of a possible 
W condensate. 

Equation (|3^) is based on massless perturbation theory. We know, however, 
that the symmetric phase is formed by massive bound states instead. This 
information may be taken into account in fitting our potential data. We shall 
also consider the ansatz (K is a modified Bessel function) 

V(m, R) = V 2 - ^-K (mR) + c 2 R (33) 

with some additional fit parameter m representing the mass of the exchanged 
particle. It should come out close to the measured W mass. 

ri i i i i i i i i n 
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Figure 18: Example for static potential in comparison to perturbation theory 

As an example we present in FigJlN the potential for {3h — 0.3434, (3q = 12 
and Mtf = 70 GeV from a run on a 3Cr lattice. The data for V(R) are obtained 
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from exponential fits to the Wilson loops 



W(R/a,T/a) = C(R)exp(-V(R)T) , 2a < R < T - 3a. (34) 

This potential is compared to the above described ansatze using massless (to = 
0) and massive (m > 0) perturbation theory. For the second case the mass 
parameter is chosen to minimize the \ 2 value of the least square fit. For this 
particular example we obtain the parameter values Vi = 0.267, c\ = er/gf = 
0.0791^1 for the massless and V 2 = 0.339, c 2 = a j g\ = O.llSgf, m = 0.93gf 
for the massive case { l m%{g%)/gt — 0.00022). In general, the fits using massive 
perturbation theory ansatz seem to describe the data somewhat better at the 
larger distances studied. 



7.3 Temperature and Higgs mass dependence of the "string 
tension" 



In Fig. [T^ we present the string tension obtained from the ansatze for the po- 
tential in (p2) and (p3), respectively, at different to§ values. In this analysis 
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Figure 19: String tension a j ' g\ from (|3^) and ( pjj ) vs. m^(gj)/g^ 



data from the higher and smaller Higgs mass Mfj cases as well as for different 
0G are used. Most of the data are from 30 3 lattices. The closest to T c and the 
/3q = 16 string tensions are obtained from fits to Wilson loops on 48 3 lattices 
(compare Fig. 12) up to distances R = 18a. The scattering of the data points 
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is an estimate for the error of the string tension. The errors obtained from the 
least square fit are obviously too small to be reliable. 

We observe that there is no significant dependence of a j g\ on A3 [Mjj] in 
both fits. At larger mj (temperatures), i.e. deeper in the symmetric phase, 
the result of the string tension based on the fit to massive perturbation theory 
( [33l ) is (not unexpectedly) close to the value reported by Teper for pure 3c? 
SU (2) gauge theory. This lends additional support to this parametrization of 
the potential. 

For the lattice distances we could explore the expected screening behaviour 
of the static potential has not yet been observed. Still larger distances are 
difficult to study due to the numerical smallness of the Wilson loops that are 
needed. 



8 Summary 

The numerical results of this study provide evidence for the first order nature of 
the thermal phase transition in the S'f7(2)-Higgs system with Higgs masses up 
to 70 GeV. The answers will soon converge concerning the upper critical Higgs 
mass above which the thermal transition is no longer first order. 

One of the motivations for the interest in non-perturbative investigations 
of the thermal electroweak phase transition in general, was to be able to check 
the reliability of perturbative calculations of the effective potential. The corre- 
sponding results of two-loop perturbation theory are confirmed for the masses 
and the renormalized Higgs condensate in the broken phase. It turns out, that 
the effect of wave function renormalization cannot be ignored. Physics that 
depends only on the potential in the vicinity of the broken minimum, can be 
systematically improved by higher order perturbation theory. 

The surface tension which is sensitive to the barrier shape of the effective 
potential for Higgs field below the order parameter at T c is systematically over- 
estimated in the perturbative calculation available so far. But this is not ex- 
clusively a problem for perturbation theory. The surface tension is notoriously 
hard to measure for weak transitions as at Mh — 70 GeV. The equilibrium con- 
figurations can be inspected and show rather structured interphase surfaces. In 
view of this, estimators of the surface tension like ours seem to be not extremely 
justified either. 

To check the viability of the dimensional reduction program in the vicinity 
of the electroweak phase transition has been another motivation for 3d lattice 
studies. The reduced model does not only make very precise predictions, but 
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seems to be reliable in the range of Higgs masses which was of particular interest 
in this investigation. In order to learn about the necessity to include higher 
dimensional operators into the effective lattice action one should explore the case 
of smaller or very much heavier Higgs masses, but this is much less interesting, 
phenomenologically. In any case, 4c? anisotropic lattices will be important for 
obtaining results (and to study masses) near to the continuum limit in physically 
large 3-volumes. 

Phenomenologically, the Standard Model is ruled out as an arena for baryon 
asymmetry generation at the electroweak transition. The 3d lattice approach 
promises to be applicable as an effective formulation of nonstandard extensions 
as well. The phase structure of the model, once fully reveiled by 3d lattice 
simulations, as well as the quantitative characterization of the strength of the 
phase transition, will be useful to inquire a multitude of 4c? extended theories 
by Monte Carlo parameter exploration. For this purpose, it is now possible to 
go beyond the easiest perturbative formulae. 

We have turned the equal weight criterion in conjunction with the multihis- 
togram Ferrenberg-Swendsen interpolation technique into a valuable iterative 
procedure. In our variant to apply this criterion, we do not need to choose any a 
priori known cut between the phases nor to associate individual configurations 
uniquely to one of the pure phases or to the two-phase mixed configurations. 
Finally, the method yields the thermodynamical weight of the latter ones. For 
any lattice size where the pure phase histograms do not overlap, the critical 
hopping parameter is obtained within the cone of infinite volume extrapolations 
of all other, more standard criteria to find /3h c - Moreover, the iterative search 
for (3hc provides the thermal energy of both pure phases in the vicinity of the 
transition, allowing to obtain an independent estimator for the latent heat. 

We have put much more emphasis than before to the properties of the sym- 
metric phase. Within the 3c? approach, information on the spectrum of particle 
like states can only be accessed through the 3c! correlation lengths. The pat- 
tern obtained lacks deeper understanding. Despite some technical progress in 
resolution of the spectrum of the 3c? transfer matrix, the interplay between the 
confining properties of the 3c? effective theory at high physical temperature and 
the space-time structure of physical excitations needs further investigations, as 
well as the nature of this confinement itself for the 3c? pure gauge theory and in 
the presence of scalar matter fields. 
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